Crossed conductance in FSF double junctions: 
role of out-of-equilibrium populations 
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We discuss a model of Ferromagnet / Superconductor / Ferromagnet (FSF) double junction in 
which the quasiparticles are not in equilibrium with the condensate in a region of the superconductor 
containing the two FS contacts. The role of geometry is discussed, as well as the role of a small 
residual density of states within the superconducting gap, that allows a sequential tunneling crossed 
current. With elastic quasiparticle transport and the geometry with lateral contacts, the crossed 
conductances in the sequential tunneling channel are almost equal in the normal and superconducting 
phases, if the distance between the FS interfaces is sufficiently small. The sequential tunneling 
and spatially separated processes (the so-called crossed Andreev reflection and elastic cotunneling 
processes) lead to different signs of the crossed current in the antiparallel alignment for tunnel 
interfaces. 

PACS numbers: 74.50.+r,72.25.-b 



I. INTRODUCTION 

Transport properties of multiterminal hybrid struc- 
tures involving a superconductor (S), connected to sev- 
eral ferromagnets (F) or normal metals (N)i£ has fo- 
cused a considerable interest recently. A superconduc- 
tor is a condensate of Cooper pairs with an energy gap 
A to the first quasiparticle excitations 3 ". In FSF double 
tunnel junctions, interesting phenomena come into play 
when out-of-equilibrium spin populations can be gener- 
ated in the superconductor^* 5 **?* 7 ^. For instance out-of- 
equilibrium effects have a strong influence on the value 
of the self-consistent superconducting gap of a FSF tri- 
laye r 4 i 5 i 6 i 7 i 8 , that can be controlled by an applied voltage. 

Cooper pairs, being bound states of two electrons 
with opposite spins, have a spatial extent £ given 
by the BCS coherence length. The limit of equilib- 
rium transport in FSF double junctions where the dis- 
tance R between the contacts becomes smaller than 
^9,10,11,12,13,14,15,16,17,18,19,20,2^ has been intensively in- 

vestigated recently in connection with the determination 
of the so-called "crossed conductance" . A voltage-biased 
crossed conductance experiment similar to the ones by 
Beckmann et aZ. 20 in the geometry with lateral contacts 
on Fig. ^ consists in measuring the current I a in elec- 
trode "a" in response to a voltage Vb — V on electrode 
"b" while a voltage Vs is applied on the superconductor. 
The trilayer geometry with extended interfaces and with 
tunnel contacts was used in a recent experiment by Russo 
et almli The crossed conductance is defined bjiii 



dl 

Ga, b (V a> V b ) = —±(V a ,V b ), 
ClVb 



(1) 



and we focus here on the case V a =Vs. Since one voltage 
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can be chosen as a reference we use V a = Vs = 0. The 
proposed interpretation of the experiment by Beckmann 
et al~* and Russo et al~ involves crossed Andreev re- 
flection and elastic cotunneling 9 ! 10 . 11 ! 12 ! 13 ! 14 . 15 ! 16 ! 17 ! 18 ! 19 , 
corresponding to transmission over two spatially sepa- 
rated contacts, in the electron-hole and electron-electron 
channels respectively, without out-of-equilibrium spin 
populations. 

Beckmann et alm& already explained that the mag- 
netoresitive effects at temperatures comparable to the 
transition temperature of the superconductor could be 
explained by charge and spin imbalance* 2 ^ due to out-of- 
equilibrium spin populations in the superconductor. This 
suggests two different explanations of the magnetoresis- 
tive effect: one explanation close to the superconduct- 
ing transition temperature based on out-of-equilibrium 
populations, and one explanation in the superconducting 
state, based on spatially separated processes at equilib- 
rium (the so-called crossed Andreev reflection and elastic 
cotunneling processes) . The goal of our article is to start 
from a description of magnetoresistive effects in the nor- 
mal state, based on out-of-equilibrium populations, and 
include superconducting correlations. We show that out- 
of-equilibrium populations can even play a role in subgap 
transport in the geometry with lateral contacts. 

More specifically we assume the existence of a region S' 
in the superconductor where the quasiparticles are not in 
equilibrium with the condensate (see Fig. For highly 
transparent FS interfaces, the geometry is characterized 
by a parameter r = Np/Ns, where Ns and Np are the 
number of channels involved in the SS' and SF contacts 
(see Fig. [TJ. The superconductor is characterized by 
a residual density of states within the superconducting 
gap pNT]/A, where pn is the normal density of states, 
and r\j A is a phenomenological dimensionless parameter 
estimated around 10~ 2 in Ref. 0> an d around 10 -4 in 
Ref. |2f| This small density of states within the super- 
conducting gap allows a new conduction channel (sequen- 
tial tunneling) that is complementary to the channels of 
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FIG. 1: Schematic representation of a FSF double junction 
with lateral contacts, with an out-of-equilibrium region S\ A 
voltage V a = is applied on the ferromagnet "a", a voltage 
Vb = V is applied on the ferromagnet "b" , and voltage Vs = 
is applied on the superconductor S. We note Ns the number of 
channels between S and S', and Nf the number of channels 
in each contact between the superconductor and the ferro- 
magnets. Different length scales are indicated on the figure: 
the spin-flip length l s f the superconducting coherence length 
£ and the energy relaxation length Ie- Other situations, such 
as the one with Ie being the smallest length scale, are con- 
sidered in the text. The situation on the figure corresponds 
to a distance R between the contacts such that f < R < l s f, 
in which case there is a finite crossed current but no magne- 
toresistive effect. The size of the out-of-equilibrium region is 
discussed in Section fff Al 



spatially separated processes. We show that the voltage 
dependence of the sequential tunneling crossed conduc- 
tance depends qualitatively on whether 77/ A C r < 1 or 
r <C 77/ A <C 1. Given the orders of magnitude of 77/ A2i2& 
and the geometry of the experiment by Beckmann et 
almQ, we conclude that the relevant regime is r <C 77/ A. 
In this case, the sequential tunneling linear crossed con- 
ductance in the superconducting state is almost equal to 
the linear crossed conductance in the normal state, which 
is apparently compatible with experiments^. 

The article is organized as follows. Preliminaries arc 
given in section ^ The cases of normal and supercon- 
ducting states with tunnel junctions are presented in sec- 
tion IIIII Numerical simulations with weak and strong 
energy relaxations and arbitrary interface transparencies 
are presented in section IIVI Concluding remarks are 
given in section [V] Some details are given in the Ap- 
pendix. 



with the two ferromagnets "a" and "b" (see Fig. ^) . The 
out-of-equilibrium quasiparticle populations in the su- 
perconductor are expected to decay over a length scale 
Iq. The length Iq can have an origin intrinsic to the su- 
perconductor, in which case it can be identified to the 
smallest between i) the spin- flip length l s f , ii) the inelas- 
tic scattering length Z e „ e , iii) the recombination length 
Ir after which quasiparticles recombine to form Cooper 
pairs^S. The length scale lo can also be due to the in- 
verse proximity effect with the ferromagnets. The den- 
sity of states induced in the superconductor due to the 
inverse proximity effect at a FS tunnel interface is propor- 
tional to 1/R 2 exp (— R/£) for a ballistic system, and to 
1/iiexp (—R/£) for a diffusive systemic, where R is the 
distance to the FS contact and £ the coherence length. 
There can thus exist a small but finite density of states 
decaying algebraically up to distances comparable to the 
superconducting coherence length. Treating the full spa- 
tial dependence of the quasiparticle populations and den- 
sity of state goes beyond the scope of our article. Instead 
we replace them by a step- function variation so that there 
exists a region S' with uniform non vanishingly small 
quasiparticle potentials, connected to the remaining of 
the superconductor S in which the quasiparticles are in 
equilibrium with the condensate. 

The number of channels Ng connecting the out-of- 
equilibrium region S' to the remaining of the supercon- 
ductor S is supposed to be large enough so that the phase 
and the chemical potential of the condensate are identi- 
cal in S and S', but the quasiparticle populations can 
be different in S and S'. Moreover the SS' contacts are 
highly transparent. Each of the FS interfaces is supposed 
to contain Np channels. The SF contacts can be highly 
transparent, like in the experiment by Beckmann et aim*, 
or have a small transparency. 

The relevant parameter characterizing the geometry of 
the FSF double junction is r = Np/Ns- The trilayer ge- 
ometry used by Russo et alS^ corresponds to r > 1 while 
the geometry with lateral contacts used by Beckmann et 
alm& corresponds to r <C 1. 



B. Hamiltonians 



The superconductor is described by the BCS Hamilto- 



HbCS = ~* (4 )0 - C /3,ff + C p,<T C «,<r) ( 2 ) 

(a,P),o 

+ A J2 ( c tl c tl + C «,I C «,T 



II. PRELIMINARIES 
A. Out-of-equilibrium region 

We suppose the existence of an out-of-equilibrium re- 
gion S' in the superconductor that contains the contacts 



The ferromagnets are described by the Stoner model with 
an exchange field h cx : 

Wstoner = ^ ~* {^a^P^ + Cp^Ca,*) (3) 
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where we supposed for simplicity that the bulk hopping 
amplitudes t are identical in the superconductor and fer- 
romagnets. We introduced in Eqs. J5J and J3J a cubic 
lattice with discrete "sites" labeled by a and f3. The 
symbol (a, (3) denotes neighboring sites on this cubic lat- 
tice, and a \ is the component of the spin along the 
z axis. The contact between the superconductor and fer- 
romagnet "a" is described by a tunnel Hamiltonian with 
a hopping amplitude t a : 



n> = t„ 



k 



"a k + c a,k Ca k , 



(4) 



where the summation runs over all sites at the interface. 
A site afc on the superconducting side of the interface 
corresponds to a site au in the ferromagnet "a" . An ex- 
pression similar to Eq. (0} is used at the interface with 
the ferromagnet "b" . The interface transparencies are 
parameterized by T a ^ = nt at bPo, where po is the normal 
density of state Pat if the superconductor, taken equal to 
the density of state pf of the ferromagnet without spin 
polarization. Highly transparent interfaces correspond to 
T a,b = 1- The transparency of the SS' contact is such that 
rs = 1. 



C. Green's function method 

1. Transport formula 

The currents are obtained by evaluating the advanced 
(G A ), retarded {G R ) and Keldysh (G+'") Green's func- 
tions2i24. The Green's functions in the superconduc- 
tor and ferromagnets correspond to the continuum limit 
since we are interested only in energies close to the Fermi 
energy. Nevertheless we introduce a discrete lattice to 
define the tunneling term Q of the Hamiltonian. 

The advanced and retarded Nambu Green's functions 
are obtained by inverting the Dyson equation, that, in a 
compact notation, takes the form G = g+g^)'S<S>G, where 
g is the Nambu Green's function of an isolated electrode, 
£ is the self-energy corresponding to the tunnel Hamilto- 
nian (0} and G is the Green's function of the connected 
structure. The symbol ® denotes a summation over the 
lattice sites involved in the self-energy. Transport prop- 
erties are obtained by evaluating the Keldysh Green's 
function G + ~ given by the Dyson-Keldysh equation 



G 



[I + G R (g>£}(g>g A 



[I- 



G 1 



(5) 



The current in the sector S z = 1/2 (corresponding to 
a spin-up electron or a spin-down hole) flowing between 
the two lattice sites a and a, is given by 

Ia, a = ^J diuTr { [i a , a &+£(u>) - t a , a Gt;-(u>)\ a z } , 

(6) 

where a z is one of the Pauli matrices, and the trace is a 
summation over the "11" and "22" Nambu components, 



corresponding to spin-up electrons and spin-down holes 
respectively. 

The Andreev current between the ferromagnet "a" and 
S' is vanishingly small since the voltage eV a is equal to the 
pair chemical potential in S'. The quasiparticlc current 
is finite because of the non equilibrium populations in 
S'. The current due to the spatially separated processes 
between S and the ferromagnet "a" is vanishingly small 
since S and the ferromagnet "a" are in equilibrium, with 
the same chemical potentials. 



2. Local Green's functions 

The local advanced Green's function of a ferromagnet 
with a polarization P is given by 



g = inp F 



1 + P 
1 - P 



(7) 



We discard the energy dependence of the ferromagnet 
Green's functions since we consider energies much smaller 
than the exchange field. 

The local advanced Nambu Green's function gs of an 
isolated superconductor takes the form 



y/A 2 - (u> - irif 



-U> + 17] 

A 



A 

-oj + irj 



, (8) 



where lo is the energy with respect to the equilibrium 
chemical potential, pn is the normal state density of 
states, and rj is a small phenomenological energy relax- 
ation paramete»2^i2£i2& ; the origin of which can be in- 
trinsic to the superconductor (inelastic electron-electron 
interaction that dominate over inelastic phonon processes 
at low temperature^) or extrinsic (the inverse proximity 
effect). We have r\ — r) cxt + ?7int, with ?y; n t the intrinsic 
value of rj, and ryext the extrinsic value. The parameter rj 
estimated to 77/ A = 10~ 4 was introduced recently^ as a 
limitation to the cooling power of microfridges based on 
NS junctions. The estimate rj/ A — 10~ 2 can be found 
in Ref.0 We will use 77/A = 1CT 2 and 77/A = 1(T 3 in 
what follows. The final results do not depend crucially 
on the precise of the absolute value of 77/A, but rather 
on how rj/A compares to r. The density of states at zero 
energy is 



Ps(uj = 0) = -lm[g 

7T 



1,1 



(w = 0)] ~ p NV /A, (9) 



reduced by a factor 77/A compared to the normal state 
density of states. 



III. CROSSED CURRENT IN THE NORMAL 
AND SUPERCONDUCTING STATES 

We start by describing magnetoresistive effects in the 
situations of linear response where either a voltage eVb <C 
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A is applied on the ferromagnet "b", or S and S' are in 
the normal state. We note t f the energy relaxation time, 
Td the transport dwell time (the average time spent by 
a quasiparticlc in S') and r s f the spin-flip time. We ex- 
amine the two cases t f <C Td <C r s f and Td <C te <C t s /, 
as well as the case of strong spin- flip r s f <C Td , te • Even 
though not directly relevant to the experiments by Beck- 
mann et a/—, we examine also the hypothesis te <C Td 
that was used recently in the study of spin imbalance 
in the FSF trilayer 4 i 5 i 6 i 7 i 8 . The transport dwell time r d 
is larger for small interface transparencies in the trilayer 
geometry, so that Td can possibly exceed te in this situ- 
ation. 

In the normal state the spin-tr current from the ferro- 
magnet "a" to S', the ferromagnet "b" to S', and from S 
to S' are given by 



from what we deduce 



I 



a~*S' 

r (<0 



Is->s> 



-NpT^ J du \fP(u) - ti f {lo - eV a ) 



^N F T^ / (Lc 



(°0 



f s /(u>) -n F {u)) 



(10) 



where f S i{u>) is the spin-cr distribution function in S', 
n F {uS) is the Fermi distribution function at zero temper- 
ature, and where the transmission coefficients Ta a \ 
and Ts are supposed to energy-independent (the full en- 
ergy dependence will be treated by numerical simulations 
in section [Wj. 



A. Weak energy relaxation (r<j < < r 3 f) 

Assuming spin conserving elastic incoherent transport 
(corresponding to Td <C te <C T sf), we impose current 
conservation for each energy to obtain 



■n F (u> - eV a ) 



(11) 



n F (u — eVb) + 



T S N S 



n F (w), 



-~n 2 k 2 F l e^exp (2R/0, 



(14) 



where £ is the superconducting coherence length, and 
where we supposed ?7i n t -C ?7oxt We obtain a cross-over 
from r <C for R <c TVs (corresponding to a point in 
the superconductor close to the contacts) to r ~ at the 
cross-over R ~ N$- The exponential increase for R > £ 
is cut-off by the intrinsic value ryj n t/A of rj/A, not taken 
into account in Eq. I|14|). 

The total currents flowing from a to S' in the parallel 
(P) and antiparallel (AP) alignments are given by 



T tot,P 



2tN F + T S N S 



T tot,AP 



2TN F + T S N S 

n 

2Tt 



(T + t)N F + T S N S 



(15) 



Vb, (16) 



where T and t denote the transmission coefficient of ma- 
jority and minority spins. We supposed r <C in the 
derivation of Eqs. lfTS|l and (|To|) . The crossed current 
is negative, and larger in absolute value in the parallel 
alignment. 

The transmission coefficients T and t are both propor- 
tional to 77/ A, while Ts is proportional to [ri/ A)(rj- mt / A). 
Assuming that ?7i n t and ?7 ex t have roughly the same or- 
der of magnitude, we conclude that the factors of order 
77/ A simplify between the numerator and denominator 
of Eqs. (H5 and JUJ in the limit T S N S > TN F . As 
a consequence, the crossed current takes approximately 
the same value in the situations where the superconduc- 
tor is in the normal and superconducting states, which is 
compatible with the experiments in Ref. l2fl 



B. 



Strong energy relaxation (te -C Td <C T s f) 



with D( CT ) 



(Ta\ a) +T^>)N F +T s Ns. 
Considering a geometry with lateral contacts used by 
Beckmann et aim*, and the magnitude of t^/Am^, Ns 
is so huge that T$Ns 3> TN F , both in the normal and 
superconducting states. This means r = N F /Ns <C r*, 
where the cross-over value of r is given by 



t 2 s Vint 



Vext 



t 2 
L F 



(12) 



Taking an estimate of ry oxt /A in the tunnel limit for 
the FS contact and in the diffusive limit for the 
superconductor—, we obtain 



A 



N F 



t F PFPN 



exp(-2iZ/0, 



(13) 



If we suppose strong energy relaxation (t f <C T d <C 
T s f), the quasiparticle distribution functions in S' are 
given by the Fermi distribution with quasiparticle po- 
tentials Vj„ and V* p for spin-up and spin-down electrons. 
The total current can then be calculated as an integral 
over energy of the transmission coefficient, and Kirchoff 
laws can be imposed on the integrated current. The spin- 
cr quasiparticle potential V° is given by 



V. 



T?N F 



qp 



{T° a + T£)N F + T S N S 



V b , 



(17) 



and the current flowing from S' to the ferromagnet "a" 
in the parallel and antiparallel alignments are given by 
the same expressions as in section Till Al 
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x ab =0.1, r/r.=10 T a b =0.1, r/r.=10" 1 x a b =0.1, r/r.=10 




FIG. 2: Variation of the normalized crossed current I a (14)/ A as a function of the normalized voltage eVb / A applied on electrode 
"b", for P = 0.5. The solid line (red) corresponds to the antiparallel alignment and the dashed line (green) to the parallel 
alignment. The absolute value of the crossed current is larger in the parallel alignment. Panels (a), (b) and (c) correspond 
to r a ,b = 0.1, and r/r, = 10 3 (a), r/r, = 10 (b), and r/r* = 10 _1 (c). Panels (d), (e) and (f) correspond to r a< b = 1, and 
r/r, = 10 3 (d), r/r, = 10 (e), and r/r, = W' 1 (f). We use rj/A = 10" 3 . 



C. Strong spin flip {r sf < r E ,r d ) 

Increasing the distance between the ferromagnetic elec- 
trodes tends to increase the transport dwell time, that 
can become larger than the spin flip length. The magne- 
toresistive effect in the crossed current decays exponen- 
tially as a function of the distance between the contacts, 
on a length scale set by the spin-flip length. The spin- 
flip length is reduced by superconducting correlations 28 , 
so that the crossed current in the superconducting state 
is reduced, as compared to the normal case. This effect 
is compatible with experiments^. 

In the case of strong spin-flip with energy conserva- 
tion (r s f «C Td <C Te), and without energy conservation 
{r s f <C te <C Td), the crossed current takes the form 



V 



-2—Nl 



T n T h 



{T a +T b )N F + T s N s ' 



(18) 



with T a = (Tj + Ti)/2 and T b = (T fc T + T b L )/2. The 
crossed current due to out-of-equilibrium populations is 
negative, but remains finite, even though there is no mag- 
netoresistive effect. The crossed current due to the spa- 
tially separated process tends to zero in the limit where 
the distance between the contacts is large compared to 



the superconducting coherence length. This may be used 
in experiments to determine whether the large distance 
behavior is due to out-of-equilibrium spin populations, 
or to the spatially separated processes. 



IV. NUMERICAL RESULTS 



A. FSF double junction with weak energy 
relaxation 



In the case of weak energy relaxation (Ta *C te <C t s /) 
we describe S' by a distribution function fgi(uj), sup- 
posed to be uniform in space within S', and determined 
in such a way as to impose current conservation for each 
energy, similarly to section UlI Al Within the numerical 
approach we can treat the full voltage dependence of the 
transmission coefficients, for arbitrary interface trans- 
parencies. The quasiparticle transmission coefficients de- 
duced from Ref. |2J are given in Appendix lAl 

The variations of the sequential tunneling crossed cur- 
rent IaiVh) in the parallel and antiparallel alignments are 
shown on Fig.|21for r a j = 0.1 and r a .h = 1. The absolute 
value of the crossed current is larger in the parallel align- 
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FIG. 3: (a) Variation of the normalized quasiparticle potentials Vq P (Vb) /A and Vq P (V b ) /A as a function of Vb/A, for r]/A — 10 
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correspond to Vq P {V b ) in the parallel (P) alignment (red, □), and to Vq P (Vi 



2 (b), and r/r, = 10 4 (c). The different curves 
in the P alignment (green, ■). V q \{Vb) is larger 



than Vq P {Vb). On panel (a) we have shown to Vq P (Vb) in the antiparallel (AP) alignment (blue, o), and to Vq P (Vb) in the 
antiparallel (AP) alignment (violet, •). We shown on the P case on panels (b) and (c) since Vq P P (V b )/A ~ Vq P AP (V b ) / A and 
Vqp AP (Vb)/A ~ V q l f (V b )/A in these cases. 
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FIG. 4: (a) Variation of the crossed current I a (V b ) as a function of Vb/A, for Tj/A = 10~ 2 , r a>b = 0.1. The three panels 
correspond to r = 1 (a), r = 10 -2 (b), and r — 10 -4 (c). The different curves correspond to I a (V b ) in the antiparallel alignment 
(red, □), and to I a (Vb) in the parallel alignment (green, ■). The absolute value of the crossed current is larger in the parallel 
alignment. 



merit than in the antiparallel alignment (which coincides 
with the normal state behavior). For small values of r/r, 
the crossed conductance for eVb > A is almost equal to 
the crossed conductance for eV& < A, in agreement with 
the argument given in section ITlI Al 



B. FSF double junction with strong energy 
relaxation 

The case of strong energy relaxation (te <C tj, -C t s j) 
can be obtained with low transparency interfaces in the 
trilayer geometry since the transport dwell time can be 
sufficiently long in this case. This case is likely to be 
irrelevant to the experiments by Beckmann et alm^ in- 
volving highly transparent FS interfaces, but is of inter- 
est for future experiments with tunnel interfaces in the 



trilayer geometry. The quasiparticle potentials Vj p and 
V~£ are calculated self-consistently in such a way that the 
integrated current satisfies Kirchoff law separately in the 
spin-up and spin-down channels. For the sake of gener- 
ality we do not treat only the case r/r, > 1, relevant to 
the trilayer geometry, but use also r/r, <C 1. 

The variation of the normalized self-consistent quasi- 
particle potentials eVJ p /A and eV} p /A as a function of 
V/ A are shown on Fig. [21 for different values of r and for 
small interface transparencies {r a ,b — 0.1). For r/r, = 
1, the normalized quasiparticle potentials eV^ p /A and 
eVjp/A increase from to a value close to unity as eV& 
increases from to A. The quasiparticle potentials are 
much reduced as r/r, decreases, due to the fact that the 
superconductor S tends to reduce the quasiparticle po- 
tentials in S' as the number of channels of the SS' contact 
increases. The quasiparticle potentials for r/r, = 10~ 2 
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and rfr-jf = 10 4 are mostly determined by the current 
flowing from S to S' and from S' to the ferromagnet 
"b". We thus obtain V^ p (V b )/A ~ V^ AP (V b ) / A and 
Vqp AP (Vb)/A ~ V q \; p (V b )/A, where P and AP corre- 
spond to the parallel and antiparallel alignments. These 
relations are due to the fact that the orientation of the 
ferromagnet "b" is reversed when going from the parallel 
(P) to the antiparallel (AP) spin orientation. 

The variation of the crossed current I a {Vb) as a func- 
tion of Vb is shown on Fig.0]for the same parameters as on 
Fig. 13 The subgap crossed current becomes very small 
for small values of r/r*, due to the fact that the self- 
consistent quasiparticle potentials are also very small. 
The absolute value of the crossed current is larger in the 
parallel alignment than in the parallel alignment. 



V. CONCLUSION 

To conclude we have discussed a model of sequen- 
tial tunneling crossed conductance based on out-of- 
equilibrium spin populations in the superconductor. The 
case of strong energy relaxation is expected to be re- 
alized for small interface transparencies in the trilayer 
geometry^SiSiL^, where the transport dwell time can be 
larger than the energy relaxation time. In the other case 
of highly transparent interfaces corresponding to the ex- 
periment by Beckmann et alm&, the transport dwell time 
is expected to be much smaller so that elastic transport 
can apply. The geometrical effects are encoded in a pa- 
rameter r, very small in the geometry with lateral con- 
tacts of the experiment by Beckmann et al. 20 . Another 
small parameter is 77/ A, proportional to the residual den- 
sity of states within the superconducting gap. There ex- 
ists a cross-over value r» such that for r <C r* the crossed 
conductance is almost voltage-independent in the elas- 
tic model. The values of 77/ A used in the literatur o 24 i 26 
indicate that the condition r <C r* is verified in the exper- 
iment by Beckmann et al. 20 . It is found experimentally 
that for the smallest distance between the ferromagnets 
the crossed conductance is the same in the superconduct- 
ing and normal phases, which can be successfully repro- 
duced by the model of sequential tunneling. Therefore a 
very small residual density of state within the supercon- 
ducting gap can lead to a sequential tunneling current in 
the geometry with lateral contacts, compatible with ex- 
periments, which is our main conclusion. In the regime of 
strong energy relaxation, not relevant to the experiments 
by Beckmann et alm&, the crossed conductance above the 
superconducting gap is much larger than the crossed con- 
ductance below the superconducting gap. Contrary to 
the elastic case, there exists a peak in the crossed con- 
ductance for eVb — A (corresponding to a large slope in 
the crossed current I a (Vb) on Fig.^J. 

In the case of tunnel interfaces, the crossed conduc- 
tance due to the spatially separated processes is posi- 
tive in the antiparallel alignmenljiiiiS while the sequential 
tunneling crossed conductance is negative. Tunnel inter- 



faces would thus constitute an experimental test of the 
possible effects. Replacing the ferromagnets by normal 
metalsSi with tunnel interfaces in the geometry used by 
Beckmann et al. 20 would also be of interest. On the the- 
oretical side it would be useful to investigate the spatial 
dependence of the out-of-equilibrium phenomena, and in- 
vestigate a more microscopic model in which the size 
of the out-of-equilibrium region would be controlled by 
the inverse proximity effect. It would be also interest- 
ing to use quasi-classical theory for describing a diffusive 
superconductor—. 
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APPENDIX A: TRANSMISSION COEFFICIENTS 
OF A FS INTERFACE 

The different terms contributing to subgap and quasi- 
particle transport at a FS interface are derived in Ref . |^ 
by Keldysh Green's function methods. In this Appendix 
we just recall these results and assume non equilibrium 
distribution functions in the superconductor. 

A first term in the spin-up quasiparticle current corre- 
sponds to transmission without branch crossing: 



0= l T ^oM W{»>)-n P (u-eV a )(u>) 
with 



r(T), 



did, 

(Al) 

2 

(A2) 

A second term in the quasiparticle current corresponds 
to transmission with branch crossing 27 : 



^SlocM = ^aT&y&M |1 + t a G^ A ^)\ 



with 



i(2) 



r(T)/ 



du>, 



(A3) 



Tl^ocH = 47r 2 ^p^>) \t a G^ A {u)\< . (A4) 

A spin-up electron from the ferromagnet "a" is transmit- 
ted in the superconductor while a Cooper pair is anni- 
hilated in the superconductor therefore producing a net 
transfer of a spin-down hole in the superconductor. 
The third term in the quasiparticle current is given by 



7 qp,eT = / T e,Ta°cH i n F ( W ~ eV ^) ~ Tl F (u + eV a )} duj, 

(A5) 



8 



with 

lSto>) = "4- 2 ^^^>^(-) (A6) 
x Re{[l + < Q G^H]G^H}. 

The density of state prefactors in the first term of 
Eq. are given by /0 o ,TPa,i/°a,i if we use G l A a A = 

— iirtaP^'lG^]^. This process therefore corresponds to 



the transmission of a spin-up electron from the ferromag- 
net to the superconductor. At the same time a particle- 
hole excitation is created at the interface, the spin-down 
hole is backscattered in the ferromagnet and the spin- 
up electron is transmitted is the superconductor. This 
process in the quasiparticle channel is reminiscent of the 
Andreev reflection term. 
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